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We revisit the light or heat-induced changes in topography of initially flat sheets of solid that elongate or 
contract along patterned, in-plane director fields. For radial or azimuthal directors, negative Gaussian curvature 
is generated - so-called “anti-cones”. We show that azimuthal material displacements are required for the 
distorted state to be stretch-free and bend-minimising. The resultant shapes are smooth and aster-like and 
can become re-entrant in the azimuthal coordinate for large deformations. We show that care is needed when 
considering elastomers rather than glasses, though the former offer huge deformations. 
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Differential growth in planar sheets [1, 2] induces topo¬ 
graphical change in order to avoid or minimize stretch en¬ 
ergy. Analogous, solid-nematic sheets are particularly rich 
since light or heat-induced length changes [3-7] are easily 
and reversibly driven. Further, their director fields n(r) that 
give the principal direction of elongation or contraction can be 
written in order to obtain any incompatibility of in-plane de¬ 
formation that is then resolved by topographical change and 
out-of-plane buckling, resulting in a desired Gaussian curva¬ 
ture. 

We revisit the simplest systems generating negative Gaus¬ 
sian curvature, namely anti-cones [8, 9] (observed by Broer et 
al [10]) that are the synthetic analogs of the ruffles, or asters, 
of Ben Amar etal[ 1,2]. In particular we examine pathologies 
arising in some modelling because of restrictions unnecessar¬ 
ily placed on shell deformations, and we present here bend- 
minimising solutions that are the optimal choice in the mani¬ 
fold of zero-stretch topographical changes from localised neg¬ 
ative Gaussian curvature. 

Effectively, heated or illuminated sheets of nematic solid 
suffer changes of metric and hence changes in Gaussian curva¬ 
ture and thus of shape [11, 12]. This metric change approach 
also has recently been exploited in such problems in an am¬ 
bitious work also concerned with general shape determination 
[13]. Our analysis is equivalent, dealing with the differen¬ 
tial geometry of space curves generating these surfaces. We 
re-examine small amplitude, stretch-avoiding solutions, since 
they point to differences with recent work [14], before we 
present large amplitude, stretch-free, bend-minimising results. 
We will also discuss the role of nematic elastomers (rather 
than glasses) that have mobile directors and can respond more 
subtly to Gaussian curvature change. 

Consider concentric circles of n in a flat sheet of nematic 
glass initially occupying the equatorial plan of the sphere in 
Fig. 1(a). Let the elongation along n(r ) be by a factor of 
X > 1 on cooling or return to the dark (or, equivalently, ra¬ 
dially arrayed n subjected to heating or light exposure with 
X < 1). There is a corresponding contraction by X v in 
the two perpendicular directions, one in-plane and the other 
through the thickness. Volume changes locally by a factor of 
Det (A) = X 1 2v where X, X v , X v are the elements of the de¬ 


formation gradient tensor A in its principal frame (based upon 
n). The opto-thermal Poisson ratio v takes values typically 
in the range (5,2) in nematic glasses [6], and v = 5 in ne¬ 
matic elastomers where volume is conserved. Note that area 
would be conserved for v — 1, as assumed by Pismen [14], 
but this value is of no particular physical significance. If elas¬ 
tic stretches are to be avoided on cooling/darkness, then one 
simply requires that, in deformation to a different topography, 
the deformation gradients take principal values X and X v cor¬ 
responding to the natural opto-thermal value under those con¬ 
ditions. We thereby automatically avoid stretch energy costs. 

In cylindrical coordinates, a reference state point is R — 
(R, <1>,Z = 0) in the initially flat disc, the equatorial plane of 
the sphere in Fig. 1(a). On cooling, its image in the target state 
is r = (p, (]), z = h ((]))), where h is the elevation from the initial 
plane (note the use of lower and upper case variables). The 
target state curve is on a sphere of radius r, see the trajectory 
in Fig. 1(a), and the in-material radius is now 

r 2 — p 2 +z 2 -A r — py^l + (h/ p) 2 . (1) 

We take a form for distortion that costs no stretch energy and 
will also minimise bend energy for small amplitudes: 

h = pA sin («(])), (2) 

with integer n for closure. The amplitude A, after scaling by 
the cylindrical radius p, is A — tana; see Fig. 1(b). The angle 



FIG. 1: (Color online) (a) An initially flat disc of radius R deformed 
into an anti-cone with an in-material radius r. Director lines are cir¬ 
cular, see for instance an n(r) indicated at radius R. (b) an n = 2 
“anti-cone”. 
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a is made between an anti-cone generator at a displacement 
antinode and the equatorial plane initially taken by the flat 
disc; see Fig. 1(a). The amplitude needed to take up the extra 
length of the perimeter with respect to the radius is a global 
requirement [8], not local, that we also return to. 

An element d.v of length in the target space image of the 
reference space element /Al<f> is 


dJ=d ‘W 


dll 

- pZ+l ^ 


We now require the sheet deforms locally according to A: 

* =!L 


R 30 


r — X v R and ds = d<I>) 

From eqn. (1) we have 

p ((()) — A~ V R/ (1 -1-A 2 sin 2 - 
/z(<]>) = A~ v RAsin(m|))/(l +A 2 sin 2 «(|)) 3 . 
Returning relations (5) and (6) to eqn (3) yields 

[l +A 2 + ( n 2 — l)A 2 cos 2 «(|)] 


ds = d<|) A ,~ V R 2 


(l + A 2 sin 2 n(| 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

(7) 


Dividing Eqn. (7) by RdO, using eqn. (4), and rearranging 
gives 



FIG. 2: Variation of the image azimuthal angle (|> with its reference 
<t> for distortion to an anti-cone for n = 2, using the infinitesimal 
form (2), but for clarity with reasonably large distortion A 1+v = 1.36, 
implying an amplitude A = 1 of the aster-like structure. The trivial 
case n = 1 of the body rotation of a disc has A = 1, no distortion, but 
the same rotation (it/4, about a diameter) is taken as that suffered by 
the antinodal generators in the n = 2 illustration. 


A = tana is identifiable only at the antinodal lines. Though 
trivial, the n — 1 case is instructive: simple geometry applied 
to the n — 1 transformation gives the mapping of the azimuthal 
angle: 


tantj) = cos a tan d>. 


( 9 ) 


gj> _i_i+v (1+A 2 sin 2 «f) 

^ < 4 > [1 +A 2 + (n 2 — l)A 2 cos 2 m))] 1,/2 

The image’s azimuthal angle (f> differs from the <I> of the orig¬ 
inal target state point. Respect for the natural distortions of 
the sheet has been built in. There are no stretch/compression 
costs, but some bend energy that we return to. 

Pismen [14] restricts all distortions in his model of anti¬ 
cones to be meridional and radial as in the more straightfor¬ 
ward case of simple cones, whereas these distortions to anti¬ 
cones described above suffer azimuthal displacements too. 
The differential method employed above is in effect equiva¬ 
lent to the metric method of Pismen, his eqns. (12)—(14), but 
with differing assumptions (about <J> and 4>) with the result 
that our anti-cones, eq. (2), do not have creases. It is inter¬ 
esting to integrate relation (8) to give ^(O), see Fig. 2, to see 
for instance the azimuthal variation in this model for an n — 2 
anti-cone which repeats at <E> — (|) — 7t. It has its first node at 
<£> = (|) = 7t/2, and its first antinode at <£> = <|) = ji/ 4. At both 
these points it is clear, for symmetry reasons, that <!> = (])■ The 
deviation of <f> from <I> is the essential difference between our 
results and those of Pismen, and is the mechanism by which 
creases are here avoided. 

The case n — 1 is simply a uniform body rotation by an 
a about a diameter since all radii at (j) and <|) + n are simple 
continuations of each other and the disc remains flat. Clearly 
X— 1 for body rotations. For n> 2 the rotation by a such that 


Alternatively, explicit integration of eqn. (8) is trivial and 
yields (9) for X = 1 and useing vT+A 2 = sec a. This analysis 
of simple rotation is a motivation for adopting li = pA sin(m|)) 
for the axial distortions into an anti-cone, and underscores a 
need to have a (f)(<J>) for non trivial cases when it already arises 
for n = 1. 

The method of [8, 9] for anti-cones was to take a global 
version of Eqn. (7) by integrating over s to give the whole 
new perimeter C, the surplus of which relative to the radius 
determines the amplitude as a function of length changes, 
A (A). Since n is along the perimeter, C — X2nR. The inte¬ 
gral of the right hand side of Eqn. (7) over (f> = (0,27t) yields 
X~ v 2nI(A,n) where the integral I is eqn (6) of [8]: 


,( "- A) =iA'Vii 


n 2 A 2 cos 2 2nnu 


+ ■ 


+A 2 sm~2nmt) 2 1 +A 2 siir 27t«M 


( 10 ) 

One sees from the current analysis that a stretch-free state is 
guaranteed globally since it is built-in locally. However one 
cannot fully specify this problem locally. From the above, 
the amplitude A and the distortion A are connected through 
I(A,n) —X l+V . As A changes, so too must the amplitude A of 
the ruffles in order that all surplus length around a perimeter 
is accommodated. The negative Gaussian curvature localized 
at the apex of the anti-cone is 2 tc( 1 — I). Eventual re-entrance 
and the transition to higher n anti-cones are discussed in [8, 9] 
and below. 
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Finite deformations: Denote by y(s) a position on a curve on 
the deformed surface with unit distance r — 1 from the ori¬ 
gin; see Fig. 1(a). It will have evolved from the circle with 
R — rX v — X v in the initial flat disc since these radii are per¬ 
pendicular to n and contract by a factor of X v . The initial 
perimeter is therefore 27tA, v which then changes by a factor of 
X (since it is along n) to 27tX 1+v . Stretch is avoided and it 
remains to minimize bend energy. 

In the Darboux frame of the curve y(s) with respect to the 
deformed surface one has 

T(s) — Y(s) ; IYI = 1 tangent 
■u(s) = u(y(s)) unit surface normal 

f(s) = tt(s)AT(s) tangent normal (11) 

with ' denoting d /d.v. The rate of change of the tangent gives 
the two components of bend: 

T'{s) — K n u(s) +Kg t(s) (12) 

where K n = u(s) ■ T' and K g = t{s) ■ T' are respectively the 
normal and geodesic curvatures. We need to minimize the 
integral of K 2 over the deformed surface. 

The surface is described by the scale-invariant curves 
r(s) = ry(s). By construction y(s) -y(j) = 1 whence ”/(s) ■ 
y(.v) -f — 0. For clarity later, we have written y = r, a unit 
vector in the surface to the curve. It is also perpendicular to Y 
and thus is f — t. From the Darboux triad we have u — rAf, 
whence from eq. ( 12 ) we have 

Kn = (rAY)y. (13) 


The total bend energy is proportional to 



where /ds is over s € (0,2nX 1+v ). The l/r with the y" arises 

because for r 1, the derivative T' is really 7 yjj = = 7 -. The 
high bend energy density at the apex of the anti-cone is dis¬ 
cussed in quantitative detail in [9], being smoothed out by 
some stretch and in any event probably not arising because 
of director escape into the third dimension during fabrication. 

The essence of the normal bend energy of the anti-cone, 

that is the \ [f Af) • factor in eq. (14) associated with 
the generator curves, can be re-written in terms of the Dar¬ 
boux frame with respect to the surface of the sphere of radius 
r — 1 on which the curve y(,v) can be also thought to live: 
T(s) = Y(s) as before, and r = u s (s) where the subscript 
s denotes “on the sphere”. The third member of the triad is 
the curve’s tangent normal on the sphere t s (s) = it s A T — 
r A Y = u. Thus the roles between the two frames have been 
interchanged: u,t -A t s ,u s . The normal curvature, eq. (13), 
on the anti-cone can be re-written as t s -T' = k), that is, it 
is the geodesic curvature of the generator curve, but on the 



FIG. 3: (Color online) (Left) The curve y(s) on the surface of a (unit) 
sphere. (Right) The ruffle surface generated by families of the the 
curves. A measure of the extra length accommodated by the ruffle 
is ratio of the (new) perimeter L to the great circle length 2 ji, here 
L/2lt = )i 1+v = 2.20, and n = 3. 

sphere (here s also denotes “on the sphere”). The normal bend 
energy of the anti-cone that we require to minimize is the same 
as minimizing f d.4 K;, 2 , which is the problem of minimizing 
the geodesic curvature of an elastica on a sphere, see Langer 
and Singer [15-18], which we now employ: 

The curve on a sphere with minimal geodesic curvature en¬ 
ergy has curvature as a function of arc length s of: 

Kg = K 0 cn(^,p) (15) 

where cn is a Jacobi elliptic function, and Ko is the maximal 
curvature along the curve. The curve must be periodic in the 
arc length, that is ^L = «4K(p) with K (p) being the com¬ 
plete elliptic integral of the first kind, since 4K(/?) is the pe¬ 
riod of cn, and L — 2nX l+v is the (new) perimeter in terms of 
the (new) radius r — 1. As before, n is the number of cycles in 
one revolution. The parameter p ensures that the curve closes 
up on the sphere’s surface - p — 0.51 in the illustration. Thus 
the maximum curvature is fixed by 

—7tA, 1+v = n4K(p). (16) 

P 

Since T' 1 = f' 2 — k | 2 + K * 2 = k | 2 + = k | 2 + 1, to plot in 

real space one needs to solve f rl = 1 + k| 2 , with k| given 
by Eq. (15) for the trajectory y(,y). Results are of the form 
of Fig. 3(left) for the trajectory y(,v) on a sphere, and of 
Fig. 3(right) for the surface generated by it. Note in particular 
that the inflection points of in k| are on a closed geodesic (here 
the equatorial great circle), there is re-entrance, there is no 
constraint of (j) = <I>, and, as expected by convexity, there are 
no sharply bent ridges where the bend energy would is large. 
For extreme elongations (large X as found in elastomers) there 
is so much extra arc length that has to be stored in the convo¬ 
luted, bent trajectory that the curves eventually intersect them¬ 
selves, forming orbit-like solutions [15-18] that are clearly 
not accessible to our deforming, impenetrable surfaces, which 
would then deviate from the solutions of Langer and Singer. 

Note that the small amplitude form adopted, Eq. (2), repre¬ 
sents the appropriate limiting case of this general solution. For 
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small amplitudes, the curves generating the anti-cone never 
deviate far from the equator of the sphere. Gauss curvature is 
a second order effect in distance and so the curve can now 
be thought of in this limit as lying instead across a cylin¬ 
der. The problem is that of periodic elastica in the plane with 
harmonic forms for displacement and curvature. To explic¬ 
itly demonstrate that the form (2) is the minimizer, first note 
that the LS minimal curvature (15) for small differences in 
length between the perimeter and the great circle, p —> 0, is 

—> Kocos (qs), taking the limit of the elliptic function for 
small modulus, and where maximal curvature also becomes 
small with p such that q = Ko/2 p is finite. The general 
form (13) of the normal curvature can be explicitly calculated 
for trajectories y(s) = (sin 0 sin(]),sin 0 cos(|),cos 0 ), in usual 
spherical coordinates, in terms of (|)(s), 0 (.v) and their deriva¬ 
tives. From the harmonic form Eq. (2), we showed 

the azimuthal variation d^/d.? is given by Eq. (7), and one 
can easily find tan0(c(>) — 1/ (A sin («()>)) and thence d0/ d<j) = 
— nAcos(nty)/ (l + A 2 sin 2 (n(j))) and thus also derivatives of 
0(s). Injecting these relations into the curvature, keeping only 
terms of order A, one obtains Kg — A(n 2 — 1) sin(/z(f)), sine 
rather than cosine appearing because of a choice of phase. 

Returning to Eq. (10) but integrating up to a (j) and an s 
rather than a complete revolution, one obtains 


A L „ 

! + T ( ” 


1 ) <>H- —s. 


Performing the complete integral, having expanded to lowest 
order in A 2 , gives an explicit form for the amplitude: 


A = 2 


/ A 1+v — 1 

n 2 — 1 


(17) 


Taking these (])(s) and A (A.) in the curvature above then gives 
the LS, minimal form = Kocos(^s), showing that the li {(.(») 
in Eq. (2) is the bend minimizing harmonic form. 

By identification, one obtains the curvature amplitude and 
the wave-vector as functions of extension A: 


K 0 = 2\J (A 1+v — 1)(« 2 — 1) q = n/ A 1+v . (18) 

One must recall that we are dealing with a unit sphere in the 
target space, which has set the scale for s and q. 

Elastomers were not discussed in [ 8 , 9] - they are more 
subtle than glasses since they can sometimes alleviate stress 
by director rotation. For instance in the azimuthal example 
considered above and by Pismen, if an elastomer disc were 
held flat on cooling so that extensile radial stress developed 
because of the deficit of natural length in the radial direction, 
then director rotation from the azimuthal towards the radial di¬ 
rection would re-attribute length from the azimuthal to radial 
direction, i.e. from a direction of (azimuthal) surplus to one 
of (radial) deficit. Little or no energy cost for such a distortion 
is required - so-called soft elasticity [19] - and the need for 
anti-cones obviated. Thus general analyses of nematic elastica 
mentioning nematic elastomers require caution. 


An elastomer example that would produce anti-cones 
would be a radial, 2-D, +1 defect being heated, that is A < 1; in 
Fig. 1(a) there would instead be radial lines of n in the initial, 
undistorted disc shown in the equatorial plane. Now a circular 
path obtains surplus length (by a factor of A ' v ) and a radius is 
in deficit (by a factor of A). The radius cannot become longer 
by rotation of the director towards it - it is already radial, and 
the glass and elastomer responses are identical in character. 
The difference is that elastomer contractions can be huge, for 
instance A —\ 0.25 on heating (—> 4 on cooling) is possible, 
and the topography changes could be accordingly larger than 
in glasses. 

We have explored and contrasted the differing assumptions 
one can make about deformations involved in anti-cone for¬ 
mation, principally the role of azimuthal displacements of 
material points. Bend energy-minimising configurations have 
been explicitly demonstrated. They embrace large deforma¬ 
tions and re-entrance. We also point to possible experiments, 
on elastomers rather than glasses, where effects are very large. 
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